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Abstract 

The’t Hooft model (two-dimensional QCD in the limit of large nnmber of colors) is 
nsed as a testgronnd for calculations of nonleptonic and semileptonic inclusive widths 
of heavy flavors based on the operator product expansion (OPE). The OPE-based 
predictions up to terms 0{l/mQ), inclusively, are confronted with the “phenomeno¬ 
logical” results, obtained by summation of all open exclusive decay channels, one by 
one, a perfect match is found. 

The issue of duality violations is discussed, the amplitude of oscillating terms is 
estimated. The method is applied to the realistic case of hadronic r decays. 


*Talk presented at the 3rd workshop on Continious Advances in QCD, University of Minnesota, 
Minneapolis, April 16-19, 1998; to appear in the Proceedings 






1 Overview 


In this talk I will present recent results obtained in collaboration with Ikaros Bigi, 
Mikhail Shifman and Nikolai Uraltsev [Q. They refer to the applications of Wilson’s 
operator product expansion (OPE) ^ to inclusive decays of heavy flavor hadrons. 
The theory of such decays is at a rather advanced stage now (see [Q and references 
therein). 

The decays of heavy flavor hadrons Hq are shaped by nonperturbative dynam¬ 
ics. While QCD at large distances is not yet solved, considerable progress has been 
achieved in this problem. The inclusive widths are expressed through OPE. The 
nonperturbative effects are then parameterized through expectation values of various 
local operators Oi built from the quark and/or gluon helds. Observable quantities, 
such as total semileptonic and nonleptonic widths of heavy hadrons Hq, are then 
given by 

i 

where q are the OPE coefficients, and p stands for a normalization point separating 
out soft contributions (which are lumped into the matrix elements {HQ\Oi{fi)\HQ)) 
from the hard ones (which belong to the coefficient functions Cj). 

There are a number of issues, both conceptual and technical, associated with 
the operator product expansion in QCD: nonperturbative contributions to the OPE 
coefficients, dependence on normalization point, violations of local duality. These 
questions are circumvented in the so-called practical version of OPE routinely 
used so far in all instances when there is need in numerical predictions. This version 
is admittedly approximate, however. The questions formulated above are legitimate, 
and they deserve to attract theorists’ attention, as they continue to cause confusion 
in the literature. 

We hnd it useful and instructive to study these subtle issues in framework of the 
’t Hooft model (see also Refs. |^). The model is 1-1-1 dimensional QCD in the limit 
Nc —^ cx). While retaining basic features of QCD - most notably quark conhnement - 
this mode is simpler without being trivial and can be solved dynamically. The theory 
is superrenormalizable, i.e. very simple in the ultraviolet domain. As a result, the 
book-keeping of OPE becomes simple, and all subtle aspects in the construction of 
the OPE can be studied in a transparent environment. 

Our study was motivated also by Ref. |]^ where heavy flavor inclusive widths were 
calculated numerically, by adding the exclusive channels one by one. It was found 
that the inclusive width THq approaches its asymptotic (partonic) value, and the 
sum over the exclusive hadronic states converges rapidly. At the same time, small 
deviations from the asymptotic value observed in the numerical analysis 0 were 
claimed to be a signal of l/mg corrections in the total width, in contradiction with 
the OPE-based result. 

We treated the very same problem analytically. The simplihcation which allows 
the analytical solution comes from taking fermions V’ (leptons or quarks) emitted in 
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Q —> g + 'ifj'il) transition to be massless. Calculating OPE coefficients and comparing 
the OPE representation (® with the sum over exclusive channels which can be found 
from the’t Hooft equation we observe a perfect match through order (1/mQ)^. 

After testing the validity of OPE, we discuss the issue of violation of local duality. 
Oscillating contributions to the total width are estimated. They are suppressed by 
the high power of l/mq which we determined. Then we apply the same method in 
the real 1+3 QCD to estimate duality violations in hadronic r decays. 

2 Setup of the problem 

2.1 ’t Hooft model 

In two-dimensional QCD the Lagrangian looks superhcially the same as in four di¬ 
mensions 

£i+i = - G%G;„ + Y.U’-V>- ™.)* . iD, = id, + Air ; (2) 

T“ denote generators of SU{Nc) in the fundamental representation and rfi the quark 
held {i is a havor index) with a mass m,; g the gauge coupling constant. 

One has to keep the following peculiarities in mind: g carries dimension of mass 
as does tpilj. With the theory being superrenormalizable no (inhnite) renormalization 
is needed; observables like the total width can be expressed in terms of the bare 
masses and bare coupling g appearing in the Lagrangian. 

The’t Hooft model is the Nc ^ oo limit of 1+1 QCD. A parameter hxed in this 
limit is ^ 

with lim = hnite . (3) 

ZTT \ NcJ Nc^oo 

This dimensionful quantity (3 which provides an intrinsic mass unit for the’t Hooft 
model, can be seen as the analog of Aqcd of four-dimensional QCD. 

The heavy hadron Hq is the bound state of a heavy quark Q and a light spec¬ 
tator antiquark The bound state is described by the light-cone wave function 
0(a;), {x G [0,1]) which is the eigenfunction of the’t Hooft equation, 

(4) 

with boundary conditions (Phq{.x = 0) = (Phq{x = 1) = 0. 

2.2 Weak interaction 

Next we need to introduce a flavor-changing weak interaction; we choose it to be of 
the current-current form: 

>Cweak = (^ 7 '"^) • ( 5 ) 




X = 


m?, - - ( 3 ^ 




+ 


X 


X 


2 








The field ip can be either the light quark or the lepton field to describe nonleptonic 
or semileptonic decays, respectively. 

For iVc —> oo factorization holds; i.e., the transition amplitude can be written as 
the product of the matrix elements of the currents q'y^Q and ip'-y^ip. For the inclusive 
widths which are discussed below the property of factorization can be expressed as 
follows: 


Mhq^Hq = j d^xi{HQ\T (6) 

= j d^xlmUf,^{x) lmT^’'{x) , 

= i{0\T {ip{x)'yf,ip{x)ip{0)j^ip{0)}\0), 

T^'^ix) = i{HQ\T{q{x)-i^Q{x)Q{d)^,q{d)}\HQ). 


This factorization follows from the fact that at > oo there is no communication 
between ipal^j.'pb and currents: any gluon exchange brings in a suppression factor 
l/Nl 

The only difference between the semileptonic and nonleptonic widths resides in 
Tl^^px), in the first case ip is the leptonic field while in the second case it is the quark 
field. At = 0, we get same (up to the overall normalization factor NP). 

Indeed, for massless fermions If^j, is given by the well-known expression. 


n^.(g) = J d^xe^^^U^,{x) = . (7) 

This expression obtained from a one-loop graph is known to be exact. For this reason 
at mp = 0 the distinction between the nonleptonic and semileptonic cases is actually 
immaterial. 

A remarkable feature of Eq. (|^) is the occurrence of the pole at = 0, which 
is specific for the vector interaction. This means that a pair of massless leptons or 
quarks produced by the vector current is equivalent to one massless boson, whose 
coupling is proportional to its momentum g^. In the case of the quark fields, it is 
known 0 from the early days of the’t Hooft model that the current ip'y^^ip produces 
from the vacuum only one massless meson, the pion. 

It means that we can use the substitution ipal^pJb ^ where the field 

(p describes a massless boson. For the Q quark decay width we have 

/^2 _ TYl^^ 

Tq = T(Q^i,i,q)=T(Q^,l,q) = - -2-i , (8) 

dTT iriQ 
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3 Operator product expansion for inclusive widths 

The l/mq expansion for inclnsive widths of heavy flavor hadrons is based on OPE 
for the weak transition operator ^ 

J d^x iT {£weak(a;)/:teak(0)} = I] O(/i)O*(p). (9) 

The local operators Oi are ordered according to their dimensions. The leading one 
is QQ with dimension dqq = 1. Higher operators have dimensions di > 1. By 
dimensional connting the corresponding coefficients are proportional to {1 / 

The coefficients q (actually, we need Imcj) are determined in perturbation theory as 
a series in {3“^/mg. These coefficients are saturated by the domain of virtual momenta 
~ mq and are infrared stable by construction. All infrared contributions reside in 
the matrix elements of the operators O^. 

At this point we see a drastic distinction between four- and two-dimensional QCD. 
In four dimensions the expansion parameter for the coefficients is the running coupling 
as{mq)] nonperturbative contributions to the coefficients coming from distances ~ 
l/mg could, in principle, show up in the form exp(—C'/ q!s(’^q)) rsj {^QCD/mqY 
where 6 is some positive index. In two-dimensional QCD such terms cannot appear: 
an analog of the exponential term above would be exp(—Cmg//?^). 

For the leading operator QQ its coefficient in the zero order in fd is Imc^g = rQ/2 
where Tg is given by Eq. (H). We found Imcgg in all orders: the only effect of high 
orders is the substitution rriq rriq — and — / 3 ^. This result is based on 

the nonrenormalization theorem we proved for the q'y^Q current at = 0. For the 
next by dimension four-fermion operators QTQqTq we show that they will appear 
only in the order. By dimensional counting it means that their contribution to 
the total widths is of order of l/rriq. 

Putting everything together we get the OPE representation for the inclusive width: 



rg 


rnp {Hq\QQ\Hq} ^ 
- 13 ^ 



( 10 ) 


Note that dependence on mq resides not only in the explicit OPE coefficient but also 
in the matrix element. The relevant matrix element is 


{Hq\QQ\Hq) 

‘^Mhq 


1 + 0 



( 11 ) 


thus the absence of linear in 1 /mq corrections is clear. We will prove a much stronger 
statement in the next section: with the accuracy of l/rug the OPE representation (pil|) 
matches the summation of exclusive widths based on the’t Hooft equation. 
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4 Match between OPE-based expressions and hadronic 
saturation 


The exclusive width T^ of two body decay Hq hn + (j) where is the mesonic [qqsp] 
state with the mass is expressed via the overlap of the initial and hnal’t Hooft 
wave functions, 


r 


n 


dvr Mhq 



dx(pn{x)(pHQ{x) 


2 


( 12 ) 


The wave function ipn(x) satisfies the same ’t Hooft equation with evident sub¬ 
stitutions, Q ^ q, Mhq Mn- 

Using this equation together with the completeness. 


Y,^n{x)^n{y) = 5{^-y) , 


we get a set of three sum rules for r„: 




d-wMHQ 

G2 


oo 


Er^ 

n=0 


AttMhq 

G2 


oo 


i: 

n=0 



^ dxcpjp^ix) = 1 , 

( 13 ) 


( 14 ) 


( 15 ) 


Note that the sums runs over all states including those unaccessible in the real 
decays of Hq, i.e. with masses Mn > Mhq- These transitions are still measurable by 
the process of inelastic lepton scattering off the Hq meson. The hrst two sum rules 
were introduced in Ref. 0. 

It is not difficult to see that the next moment, is a divergent 

sum because (p'jj (x)/x^ would no longer be integrable. It defines the asymptotics, 

W _ 

r„ oc l/M®, at large n. For this reason the “unphysical” part of the sum (|^) gives 
only l/rug corrections and we derive from the sum rule (M) 



rg 


_Mhq 



^Hq (^) + ^ 



( 16 ) 


This expression coincides with the OPE result of Eq. (ED as seen by rewriting the 
matrix element in Eq. (pUD in terms of the ground state wave function lphq{x). (The 
factor mQ/{mQ — accounts for nonlogarithmic dependence of QQ operator on 

normalization point). 

Thus, we have demonstrated the perfect match between OPE and the hadronic 
saturation for the inclusive width. 
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5 Deviations from local duality 


Having established a perfect match between the OPE prediction for the total width 
and the result of the saturation by exclusive decay modes, through Oil/mq), we 
must now turn to the issue of where the OPE-based prediction is supposed to fail. 
The failure usually goes under the name of “duality violations”, a topic under intense 
scrutiny in the current literature. 

We must precisely dehne what is meant by duality and its violations. Sum rules 
relating certain moments of the imaginary part of transition amplitude to matrix 
elements of consecutive terms in the OPE series will be referred to as global duality. 
Taken at their face value, they are exact, to the extent of our knowledge of the 
coefficient functions and matrix elements of the operators involved in OPE. Therefore, 
it does not make any sense to speak about violations of the global duality. 

The notion of local duality on the other hand requires further assumptions. The 
difference between the OPE-based smooth result and the experimental hadronic mea¬ 
surement is referred to as the duality violation meaning the violation of local duality. 
Thus, the duality violation is something we do not see in the (truncated) OPE series. 
The duality-violating terms are exponential in the Euclidean domain and oscillating 
(like sin[(E/AQCD)^]) in the Minkowski domain. 

The appearance of duality violations in the form of oscillating terms is evident in 
the’t Hooft model where the spectral density is formed by zero-width discrete states. 
Indeed, each time a new decay channel opens dTHq/dniq experiences a jump (P^q is 
continuous) , so that immediately above the threshold, dVH q !dmq is larger than the 
smooth OPE curve. In the middle between two successive thresholds it crosses the 
smooth prediction, and immediately below the next threshold dVH q/ dirtq is lower 
than the OPE-based expectation. To estimate the amplitude of oscillations we use 
the following model for exclusive widths P„ at the range of in the vicinity of Mhq '. 

, , M2 _ m2 

"y (17) 

Then for the oscillating part AP”'" of the total width P^^^ as a function of mq we get 


AP osc 




/9 

Mhq 


9 



x) 


P 

6 . 


(18) 


where 


X = Fractional Part of 




X e [0,1). 


The plot of AP°'’‘' is presented in Fig. |I|. The amplitude of oscillations is /‘2){j3 / MhqY . 

Thus, we see that violation of local duality is present but suppressed as I/uIq. If 
we could average over mq in a sufficiently large interval we would get an exponen¬ 
tial suppression. Then it would be possible to consider the OPE-based predictions 
beyond l/rug. 


6 









Figure 1: Oscillations in THq- The ratio AF^'^/G^/? is presented as a function of 


6 Duality violations in r decays 


Let us discuss along similar lines a quantity of practical interest in 1+3 dimensions, 
namely the normalized hadronic r width 


Kr 


F (r —> z/,r + hadrons) 
F(r- ^ l^rG-pe) 


ds / 

Jo M^[ 



(19) 


where p{s) is the sum of spectral densities in the vector and axial-vector channels. 

To estimate the oscillating part of R-r we consider the limit of large and M^-. 
For large Nc the spectrum of 1+3 QCD is expected to consist of an inhnite comb 

[ 0 . 


of narrow resonances - in complete analogy to the ’t Hooft model 


To keep 


the closest parallel to it we further assume that the high excitations are equally 
spaced in m?. This agrees with the general expectation of a string-like realization of 
conhnement leading to asymptotically linear Regge trajectories. The masses of the 
excited states in, say, the vector channel are then given by = m^ + 2n/a' o, with 
a' being the slope of the Regge trajectory. Experimentally one hnds 2la' ':^2 GeV^. 

Then spectral density p(s) at large values of s approaches the form: 


p{s) = 2iVc ■ ^ 5 (— - n j ; = — , (20) 

n=l V(T / a 

Equation (1^) is clearly not expected to hold at moderate and small values of s where 
the vector and axial-vector channels are drastically different and the resonances are 
not equidistant. However, details of the spectral densities at small s play no role in 
duality violation, they change only the regular terms of the 1/M^ expansion. The 
spectral density in Eq. (^O]) is dual to the parton model result; i.e., it coincides with 
it after averaging over energy, (p(s)) = 2Ac, and Rr approaches Nc asymptotically. 
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The sum over resonances in is calculated analytically: 


R-r 

^OPE 

6R°^^ 


where 


^OPB ^ ^^osc ^ (21) 

M2 ^ 30 1^M2 ) ’ 

(£)■ 


■ X{1 - X){1 - 2x) l^j + : 


x'^il — xY - 

^ ^ 30 


X = Fractional Part of 



X e [0,1) . 


The result is a sum of two functions of M^, the hrst one, is a smooth function 

expandable in 1/M2, it coincides with the OPE prediction in the model. The second 
one, 6R°^^, oscillates with the period a^, see the plot of 6R°^^/Nc in Fig. The 
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Figure 2: Oscillations in Rt-. The plot of 6R°^^/Nc is presented as a function of 
Mila‘S. 


dominant component of 6R'^^'^/Nc scales as 1/M®. This oscillating component vanishes 
at values of corresponding to the new thresholds, and at one point in the middle 
between the successive resonances; its second derivative has a jump at the thresholds, 
so we deal with a nonanalytical dependence on M^-. The average of 6R°^^ vanishes 
while the amplitude of oscillations amounts to 


_ 1 ( 

“rT „„ ^ VH \mII 


( 22 ) 


Taking our estimate of the oscillation amplitude at its face value and using the 
actual value of the r mass in Eq. (|2[) we hnd 6R°^^/Rr ~ 3%. One should not 
the number take too literally for many reasons: hrst of all the r mass is not much 














larger than the spacing between the resonances, jM'^ ~ 2/3; second, iVc is not large 
enough to warrant the zero width approximation, the nonvanishing widths lead to a 
further suppression of deviations from duality. Still we believe that the consideration 
is instructive in a qualitative aspect. 
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